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FOREWORD 


This  is  the  final  report  for  the  surface  scatter  study  program  undertaken 
by  the  Optical  Sciences  Center,  University  of  Arizona,  Tucson,  Arizona.  The 
work  was  performed  under  Air  Force  Contract  F04701-72-C-0181 .  The  program 
was  carried  out  under  the  auspices  of  the  Air  Force  Space  and  Missile  Systems 
Organization . 

The  principal  investigator  for  the  Optical  Sciences  Center  was  Roland  V. 
Shack.  The  project  engineer  and  co-investigator  was  Michael  A.  DeBell. 


SUMMARY 


Scattering  theory  is  developed  as  an  extension  of  the  angular  spectrum  of 
plane  waves  approach  to  diffraction.  Expressions  useful  in  the  evaluation 
of  scatter  are  obtained  through  the  introduction  of  random  perturbations  in 
the  pupil  function.  When  the  surface  statistics  are  known  the  expected 
value  of  the  transfer  function  can  be  found.  The  particular  case  of  a  sur¬ 
face  with  Gaussian  statistics  is  then  treated  and  expressions  are  derived 
that  predict  the  scattered  flux  distribution  for  a  surface  with  such 
statistics . 

Apparatus  is  described  that  is  capable  of  measuring  the  scattered  flux 
distribution  in  the  range  of  two  to  ninety  degrees  from  the  specular  direc¬ 
tion.  Using  this  apparatus  the  distribution  of  light  scattered  by  a  variety 
of  plane  mirror  surfaces  as  a  function  of  roughness,  angle  of  incidence,  and 
wavelength  is  measured  experimentally. 

The  intensity  of  the  scattered  light  is  shown  to  vary  with  the  inverse 
square  of  the  wavelength  when  the  total  scattered  flux  is  small  compared  to 
the  specularly  reflected  flux.  This  theory  is  consistent  with  the  data. 

A  numerical  integration  of  the  scattered  light  distribution  yields  the 
total  scattered  flux.  A  relationship  between  the  total  scattered  flux  and 
the  surface  roughness  is  derived.  The  theoretical  results  and  the  estimated 
roughness  were  found  to  be  of  the  same  order  of  magnitude. 
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INTRODUCTION 


When  light  is  incident  upon  a  surface  the  reflected  flux  can  be  described 
by  considering  two  components,  the  specular  component  and  the  diffuse  compo¬ 
nent.  The  specular  component  is  reflected  in  the  plane  of  incidence  at  an 
angle  equal  to  the  angle  of  incidence  and  is  coherent  with  the  incident 
beam.  Light  scattered  away  from  the  specular  direction  is  said  to  be 
diffusely  reflected.  Diffusely  reflected  light  is  a  diffraction  phenomena 
and  its  distribution  is  attributable  to  the  surface  properties  of  the  re¬ 
flector.  Properties  that  can  be  used  to  characterize  such  surfaces  are 
roughness,  autocorrelation  length,  and  the  statistical  distributions  of 
these  two  parameters. 

Two  special  cases  are  of  interest:  the  perfect  reflector  and  the  per¬ 
fectly  diffuse,  or  Lambertian,  reflector.  A  perfect  reflector  alters  the 
direction  and  not  the  distribution  of  the  incident  light.  A  plane  surface 
that  is  perfectly  diffuse  obeys  Lambert's  cosine  law,  1(0)  =  Iocos0,  where 
1(0)  is  the  radiant  intensity  at  an  angle  0  to  the  surface  and  ID  is  the 
radiant  intensity  measured  normal  to  the  surface.  Theories  put  forth  should 
reduce  to  these  special  case's  in  the  limit  and  it  is  in  this  regard  that  we 
mention  them. 


Previous  Theoretical  Work 


Scattering  from  surfaces  is  a  diffraction  phenomena  and  as  such  has  been 
treated  as  a  boundary  value  problem  where  Maxwell's  equations  are  applied. 
Two  categories  exist  in  this  general  body  of  theory:  that  in  which  scat¬ 
tering  is  treated  for  a  specific  model  of  the  surface  and  that  in  which  the 
surface  is  statistically  described  by  a  distribution  of  heights  and  slopes 
about  a  mean  surface. 

Examples  of  theoretical  approaches,  for  a  specific  surface  model,  have 
been  proposed  by  Beckman  and  Spizzichino  (1963),  Berreman  (1967),  Hunderi 
and  Beaglehole  (1970),  and  Twersky  (1957).  Twersky  and  Berreman  have  con¬ 
sidered  surface  models  for  which  it  is  possible  to  obtain  exact  solutions 
for  the  electromagnetic  field.  In  these  cases  the  Kirchhoff  approximations 
are  applied  to  the  boundary  conditions  in  order  to  evaluate  the  Helmholtz 
integral.  Their  models  consist  of  randomly  distributed  hemispherical 
irregularities  on  an  otherwise  smooth  surface.  Twersky  arrived  at  complete 
solutions  by  considering  only  perfect  conductors  while  Berreman  made  first- 
order  calculations  applicable  to  surfaces  of  any  conductivity.  Hunderi  and 
Beaglehole  adopted  an  approach  that  yields  the  reflectance  of  a  surface 
model  consisting  of  spheres  located  just  above  a  smooth  metal  surface. 
Although  this  model  of  a  scattering  surface  is  not  as  realistic  as  that  of 
Berreman,  it  can  be  solved  exactly  for  a  metal  with  any  conductivity. 

Theory  relating  to  periodically  rough  surfaces  is  discussed  by  Beckmann 
and  Spizzichino.  They  use  the  Kirchhoff  method  of  approximating  the  bound¬ 
ary  conditions  on  a  perfectly  conducting  periodic  surface  and  then  find  the 
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corresponding  field  by  using  the  Helmholtz  integral.  Making  additional 
approximations  they  deal  with  the  surfaces  that  are  not  perfect  conductors. 

Perhaps  the  best  known  theory  dealing  with  statistical  scatter  from  sur¬ 
faces  having  roughness  distributions  is  that  of  Davies  (1954) .  His  theory 
considers  only  the  phase  modulation  of  the  incident  light  by  the  height  var¬ 
iations  along  the  surface.  The  specular  and  scattered  light  are  found  by 
summing  the  radiation  from  Huygens  wavelets  at  the  surface.  The  height  var¬ 
iation  is  then  assumed  to  be  normally  distributed  about  its  mean  with  an  rms 
deviation  of  a.  The  scattering  theory  then  predicts  that  the  ratio  of  the 
specular  reflectivity  of  a  rough  surface  to  that  of  a  smooth  surface  varies 
as 


R/R 

o 


e- (4-rra/A) 2 


and  makes  similar  predictions  about  the  angular  distribution  of  scattered 
light.  The  approximations  made  in  this  development  limit  the  theory  so  that 
it  is  valid  only  if  the  typical  distance  separating  surface  structures  is  on 
the  order  of  a  wavelength  or  larger.  This  theory  of  Davies  was  extended  by 
Bennett  and  Porteus  (1961)  and  by  Porteus  (1963) .  Work  by  Beckmann  and 
Spizzichino  (1963)  also  extended  the  statistical  approach  of  Davies. 


Previous  Experimental  Investigations 

Lord  Rayleigh  (1901)  was  one  of  the  earliest  investigators  of  scattered 
light.  He  noted  the  effect  of  poorly  polished  surfaces  on  optical  perfor¬ 
mance.  His  work  examined  the  effects  caused  by  surface  roughness,  wave¬ 
length,  and  angle  of  incidence  on  the  reflected  beam.  The  variation  in 
brightness  with  angle  of  observation  for  various  surfaces  was  studied  by 
Jones  (1922).  He  used  what  he  termed  a  "goniophotometer";  a  device  having 
means  of  directing  a  collimated  beam  on  a  sample  surface  and  a  photometer 
arm  arranged  to  measure  the  brightness  of  the  surface.  The  angle  of  illumi¬ 
nation  could  be  varied  and  the  angle  of  viewing  held  constant  or  vice  versa. 
Slater  (1935)  improved  on  this  goniophotometer  and  was  able  to  provide  quan- 
titized  measurements  in  either  transmission  or  reflection. 

Photometric  measurements  by  Hunter  (1940)  employed  a  barrier  layer  photo¬ 
cell  to  measure  the  diffuse  reflectance  of  a  broad  range  of  samples.  His 
work  included  a  detailed  report  on  sources  of  error.  Hunter  used  aluminized 
"perfect  mirrors"  as  reference  standards.  Some  of  the  first  work  where 
scattering  itself  was  a  primary  concern  was  reported  by  Aughey  and  Baum 
(1954).  The  scattering  profiles  they  measured  led  them  to  direct  conclu¬ 
sions  about  the  physical  properties  of  the  scatterer. 

The  scattering  theory  of  Davies  (1954)  was  expanded  and  experimentally 
investigated  by  Bennett  and  Porteus  (1961).  From  this  work  and  work  that 
followed  shortly  thereafter  (Bennett,  1963  and  Porteus,  1963)  the  reflec¬ 
tance  properties  of  samples  with  a  measured  surface  roughness  were  directly 
compared  to  theory.  In  their  work  Bennett  and  Porteus  reported  good  agree¬ 
ment  between  the  measurements  on  their  samples  and  the  theory  they  had 
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developed.  The  thrust  of  these  measurements  led  to  investigations  at  the 
Optical  Sciences  Center  by  D.  McKenney,  L.  Mott,  G.  Orme  and  M.  DeBell. 

Semplak  (1965)  and  Blazey  (1967)  along  with  others  (see  annotated  bibli¬ 
ography)  used  lasers  as  a  light  source  for  their  scattering  measurements. 
Blazey,  using  a  goniometric  scanning  device  reported  results  in  good  agree¬ 
ment  with  the  Beckmann  theory.  His  measurements  were  made  on  multi-layer 
dielectric  coatings  typical  of  those  used  in  the  laser  industry. 

Orme  (1972),  using  a  clever  diffraction  technique,  was  able  to  explore 
the  zeros  of  a  single  slit  diffraction  pattern  and  evaluate  scatter  to 
within  one-third  of  a  degree  of  the  specular  beam.  Orme  made  wide  angle 
measurements  using  a  goniometric  device.  Using  this  technique  he  was  able 
to  explore  all  but  the  limited  region  of,  1.5°  to  5°,  from  the  specular  di¬ 
rection,  in  the  plane  of  incidence.  His  work,  as  had  the  work  of  Aughey  and 
Baum,  pointed  to  the  possibility  that  scatter  was  caused  by  two  different 
mechanisms.  The  inability  of  Orme  to  make  measurements  in  the  1.5°  to  5° 
range  and  the  questions  raised  concerning  multiple  scattering  phenomena  have 
provided  the  stimulation  to  undertake  the  work  reported  here. 
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II 


THEORY 


The  Angular  Spectrum  of  Plane  Waves  Approach  to  Diffraction  Theory 

The  aim  of  this  section  is  to  cast  diffraction  theory  in  a  format  that 
will  make  it  useful  in  examining  the  effects  of  scattering.  We  will  dis¬ 
tinguish  scatter  from  diffraction  by  defining  scatter  as  "diffraction  from  a 
surface  or  medium  whose  properties  are  statistical  in  nature."  Our  approach 
will  be  to  consider  a  disturbance  in  an  initial  plane,  Pj,  and  to  derive  the 
necessary  theory  that  will  allow  us  to  predict  the  field  in  a  subsequent 
plane,  P2. 

Let  the  complex  amplitude  distribution  at  the  plane  P^  be  U(x,y,0)  and 
consider  the  geometry  shown  in  Fig.  1  below.  The  scalar  disturbance  on  Pj 


will  be  considered  the  only  radiation  contributing  to  the  field  at  It 

will  be  assumed  that  the  light  propagating  from  Pj  to  P2  obeys  the  Helmholtz 
equation.  Further  we  will  assume  that  the  Fourier  transform  of  the  scalar 
field  U(x,y,0)  exists.  (We  can  make  this  assumption  since  as  Goodman  points 
out,  "physical  possibility  is  a  valid  and  sufficient  condition  for  the 
existence  of  a  transform.") 

We  can  now  define  the  Fourier  transform  relationships  that  exist  for 
planes  P^  and  P^. 
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(tz  has  a  parametric  relationship  since  it  is  a  function 
of  the  observation  plane.) 


Notice  that  the  integrand  in  Eqs.  (3)  and  (4)  has  a  form  similar  to  that  of 
a  plane  wave.  The  equation  for  a  unit  amplitude  plane  wave  is 


e 


ik  *r 


Here  k  is  the  propagation  vector  and  r  is  the  position  vector  (see  Fig.  2) . 
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where 


y 

and  z 

are 

unit  vectors ,  and 

a. 

6 

and  y 

are 

direction  cosines 

a2 

+ 

32  +  y2 

= 

1. 

(5) 

(6) 


The  unit  amplitude  plane  wave  can  then  be  rewritten  as: 


ik*  r  _  i(ax  +  gy  +  yz)2tt/A 

v  0 
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(7) 


Across  the  plane  z  =  0,  the  integrand  of  Eq.  (3)  can  be  thought  of  as  the 
intersection  of  that  plane  with  a  plane  wave  traveling  with  direction 
cosines : 


a  =  A£,  3  =  An  and  y  =  A  -  (A^)2  -  (An) 2 . 


(8) 


It  now  becomes  convenient  to  scale  the  variables  x,  y,  and  z  to  the 
wavelength,  A.  These  new  variables  are  denoted  x,  y,  and  z; 


x_ 

A  ; 


(9) 


With  this  notation,  Eqs.  (1)  -  (4)  can  be  rewritten  as: 
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(12) 


U(x,y;z) 


A(a,B;'z)  ei2’(0i  *  ^dadB. 


(13) 


Here  we  have  cast  the  integrand  into  the  form  of  plane  waves  modified  by 
Fourier  coefficients. 

The  disturbances  in  planes  Pq  and  ?2  must  also  satisfy  the  Helmholtz 
equation: 


(V2  +  k2)U  =  0. 


(14) 
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Here, 
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82 


92 
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In  the  scaled  coordinate  system,  V2  =  A2V2,  and  k2  =  A2k2  =  (2tt)2.  Wi  th 
these  substitutions,  Eq.  (14)  becomes: 

(V2  +  (2it)2)U  =  0.  (lf0 

If  we  apply  Eq.  (15)  to  Eq.  (13)  we  get: 
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hence. 
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The  Fourier  transform  of  the  function  in  brackets  is  zero  so  the  function 
itself  must  be  identically  zero.  Thus  we  can  write: 


A(a, B; z)  +  4iT2Y2A(a,B;z)  =  0, 


(19) 


since  y2  =  (l-a2-32) .  This  equation  has  a  solution  of  the  form: 


A(a,8;z)  =  Ao(a,3;0)  e 


■  i2TTyz 


(20) 


Equation  (2)  then  relates  the  Fourier  transforms  of  the  scalar  fields  in 
planes  P^  and  P2-  Indeed  Eq.  (20)  can  be  rewritten  in  terms  of  a  transfer 
function  for  free  space,  H; 


A(a,3;z)  =  Aq (a, g ;0)H(a, 8 ; z) .  H(a,g;z) 


=  e 


i2iTYZ 


(21) 


We  have  thus  far  applied  no  restrictions  on  y  and  two  regions  of  interest 
are  apparent;  that  for  real  values  of  y  and  that  for  imaginary  values. 


Y 


[1 


(a2  +  B2)]3s 


for  (a2  +  32)  <  1 


Y  is  real 


for  (a2  +  32)  >  1 


Y  is  imaginary  (22) 


Where  y  is  real  the  various  plane  wave  components  of  the  angular  spectrum 
propagate  in  different  directions  and  thus  reach  the  "observation  point"  in 
plane  P2  with  an  associated  phase  factor 


i2trYZ 

e 


When  y  is  imaginary  the  exponential 


ei2iTYZ 


becomes  a  damping  function  of  the  form  e~K'z  where  K  =  2tiV (a2+82)  -1 .  Waves 
characterized  by  these  values  of  y  are  rapidly  attenuated  in  the  z  direction. 

The  spread  function  associated  with  our  transfer  function  is  its  Fourier 
transform. 
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■  (‘-Sf) 


This  spread  function  can  be  cast  into  the  form: 


a(r)  e 


i<Kr) 


where  r2  =  x2  +  y2  +  z2 


The  righthand  side  of  Eq.  (23)  is  straightforwardly  separated  into  its  real 
and  imaginary  components,  a  +  ib; 


sin2-rrr  cos2Trr 

~~r  2frf2 


cos2irr 


sin2irf 
2irr  2 


To  represent  this  in  phasor  form  we  first  recall  the  graphical  representa¬ 
tion  shown  in  Fig.  3. 


Fig.  3.  Phasor  representation. 


a(r)  =  /a5 


XF  f1  * 


The  phase  can  be  expressed  as: 


(f>  =  arc  tan(— ) 

a. 


=  2-rrr  -  arc  tan  (2irr) 
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Hence  our  spread  function  can  be  written  as: 


a(r)  e 


i(f)(r) 


Cl  *  C2.f)2)H  eit2"  -  tan2ri) 


(28) 


The  measurements  we  wish  to  make  have  two  characteristics  that  should  be 
considered:  measurements  are  made  on  a  spherical  surface  and  they  are  made 
at  a  distance  from  the  aperture  plane  that  is  large  compared  to  the  size  of 
the  "initial  disturbance"  and  the  wavelength  of  the  light  that  is  propa¬ 
gating.  Figure  4  shows  the  geometry  that  we  wish  to  consider. 


Hemispherical 

Surface 


Fig.  4.  Relationship  of  diffraction  plane  to 
observation  surface. 


The  primed  (')  plane  is  the  pupil  plane  and  the  spherical  surface  of 
radius  ft  is  where  we  wish  to  measure  the  resultant  disturbance.  In  terms  of 
what  has  been  previously  discussed  we  can  write: 

A(a, 6 ; z)  =  Ao(a,6;0) 

U(x,y;z)  =  F{ A(a, B ; z} 
and  can  then  conclude. 


U(x,y; z) 


F(Aq  (a,  3 ;  0) 


ei2TTyz 


}. 


(29) 
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Now  by  the  convolution  theorem  we  have: 


U(x,y;z)  =  F(A  fa,B;0)}  *  F{ei2iry2} 


=  U(x' ,y' ;0)  *  (  i  - 


1  \  z  e 


i2TrR 


(30) 


2irR  /  R  R 


The  convolution  integral  can  then  be  taken  "over"  the  pupil  plane: 


U(x,y;z) 


U(x',y’;0)(i 


,  -  /  i2'frr\ 


(31) 


where 


r  =  /(x-k1)2  +  (^-y')2  +  g2. 


Now  assume  that  we  are  many  wavelengths  away  from  the  pupil  plane.  The 
term  1/2'rTr  then  becomes  insignificant  in  the  integral  compared  to  the  other 
term.  Equation  31  can  then  be  approximated  by: 


U(x,y;z)  =  i 


-  /  i2-rrr\ 

U  (x ' ,  y ' ;  0)  (|-)  (  — ~ —  )  <lx '  dy ' 


(32) 


Substituting  R  +  (r-R)  =  r,  F.q.  (36)  is  rewritten  as 


U  (x,y;  z)  =  i 


U(i.y;0)('4YI)2fS^')  ei2ir  V-R)d' 

\r)\t/  \  R  / 


’dy’.  (33) 


But 


r-R 


x’2  +  y’2  2R 


r  +  R  r  +  R 


(v’  +  &</')» 
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where 


Equation  (33)  can  now  be  expressed  in  terms  of  the  direction  cosines  and 
the  radius  R. 


u(V30;R) 


iy  e 
o 


i2rR 


u(x',9';0)(|)2  e 


i  2tt  ( — ' 2 


r  +  R 


-i2ir  — —  (a  x'  +3  y') 
r+R  0  0 

x  e  dx ' dy ' 


We  now  make  two  approximations: 


(35) 


(a)  R  +  r  =  2ft  1 


r  -  R 
2R  / 


i*  ~  R 

2R,  since  -  <<  1  and 


(b)  £ 

^  1  r 


r  +  R-r 


=  1  + 


ft-i 


(j)2  ,  !  *  (tef  =  1, 


R-r  .  .  * 

since  — <<  1 . 


(36) 


When  the  approximations  expressed  in  Eqs.  (36)  are  substituted  into  Eq.  (35) 
we  have: 


uCveo'« 


iv  e 
o 


i2irR 


R 


i  2tt 

U(x',y' ;0)e 


X 


”i2Trfa  xf 
o 
e 


+ 


V') 

dx ' dy ' . 


(37) 
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Equation  (37)  is  the  classical  Fresnel  diffraction  expression  modified  by 
a  spherical  wave.  The  further  approximation  that 


is  equal  to  1  leads  to  the  result: 


U  (aQ,  gQ  ;R) 


iy  e 
o 


i2irR 


R  J 


-i2ir(a  x'+g  y') 

U(x,,y';0)e  °  °  dx'dy'. 


(38) 


This  is  the  Fraunhofer  pattern  as  modified  by  a  spherical  wave.  From  Eq. 
(11)  we  know  that  Eq.  (38)  can  be  rewritten  as: 


i2frR 

iy  e 

U  (a  ,g  ;R)  *  — 5 -  A  (a  g  ;0)  . 

o  o’  ~  ov  o  o  1 

K 


(39) 


We  have  a  spherical  wave  modulated  by  the  Fourier  transform  of  the  distur¬ 
bance  in  the  initial  plane.  Now  let 

U(x',y')  =  a  x(x',y').  (40) 

From  Eqs.  (39)  and  (40)  we  can  write 


R  -i2rR 
- -  e 

1Yo 


U(a0,B0,R) 


A(ao’V 


J  J 

—  CO 


-i2n(a  x'+g  y') 

x(x',y')e  0  0  dx'dy'. 

(41) 


The  power  theorem  will  give  us  the  total  flux,  viz.. 


CO 


—  CO 


U(a  ,g  , R) I  2  da  dg 
o  o  1  o  o 


oo 


J 


—  OO 


x (x' ,y ') | 2  dx'dy' 


<1.  (42) 
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Noting  that  dm  =  da0dB0/Yo>  Eq.  (42)  can  be  rewritten  as: 


—  U(a  ,3  ,  R)  2  dco 
Y  1  o  o 
o 


? 

a  a 
o 

A2 


=  4>  =  total  power. 


(43) 


The  radiant  intensity,  I(w/sr),  is  then 


■(VV 


d$ 

dto 


f  lu(»o-eo’R) 
o 


2  _ 


Y  I A (a  ,3  ) I  2 . 

'o'  o  o'  1 


(44) 


We  note  that  the  radiant  intensity  is  directly  proportional  to  y  and  thus  we 
have  a  cosine  dependence- 


Scattering  Theory 


In  general  the  pupil  function  will  include  random  perturbations  that  will 
lead  to  scatter.  In  order  to  treat  these  perturbations  we  consider  the 
transfer  function  V,  which  is  the  Fourier  transform  of  the  spread  function, 
i.e.,  | A | 2  +->-  V.  The  autocorrelation  theorem  allows  us  to  write  the  trans¬ 
fer  function  as  in  Eq.  (45)  below. 


V(5' »n ') 


oo 


* 

U(xf 


j  j 
—  00 


y' 


X  U*(x' 


(45) 


oo 


—  OO 


X  T*(x'  -  ~  ,  r  -  -^~— ) dx ' dy ’ .  (46) 


If  we  assume  that  the  pupil  function  has  random  phase  variations  and  no 
amplitude  variations  (such  as  a  mirror  in  the  pupil  would  introduce)  we  can 
write  x  as: 


T 


i2ir(W  +  Wn) 

0  R' 

> 
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(47) 


where  WR  indicates  the  random  component.  Further,  if  we  introduce  the  nota¬ 
tion: 


W(+)  =  W(x  +  ,  y  +  ij— ) 

and 

W(-)  =  W(x  -  11  ,y  -  nLh 


(48) 


and  take  the  average,  we  obtain 


<V(C  ,ri  ')> 


i2ir(W  (  +  ) 
t(-)|e  ° 


woC-)) 


i 2tt  (W  (+) 
x  <e  R 


wRC-)) 

dx ' dy ' . 


(49) 


( 


The  term, 


i 2tt  (W  (+)  -  W  (-) ) 

<e  R  > 


in  Eq.  (49)  is  the  only  term  of  the  equation  that  is  random.  For  stationary 
processes,  such  as  those  which  describe  mirror  surfaces,  this  term  is  not  a 
function  of  x'  or  y'  and  may  thus  be  taken  outside  of  the  integral. 


<V(?',n') 


|t  (+)  I  |t  (-) 


i2ir  (W  (+) 

x  e 


W  (-)) 
o 


dx ' dy ' 


(50) 


The  expected  value  of  the  transfer  function  may  be  written  as: 


<■ 


i2  tt(Wr(+) 


( 

(51) 


Here  VQ  can  be  thought  of  as  a  transfer  function  for  a  surface  without 
irregularities  and  the  expected  value  of  the  term  containing  the  random  com¬ 
ponents  as  a  transfer  function  for  the  irregularities.  It  remains  to 
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determine  the  form  of  this  expected  value.  Let  us  assume  a  normal  distribu¬ 
tion  for  WR(+)  -  WR(~) .  From  random  variable  theory  we  find  that  the 
expected  value  for  two,  real  valued,  stationary,  jointly  normal,  random  pro¬ 
cesses  X(+)  and  Y (+)  is  given  by: 


where 


*(+)  +  iY(+) 


nx  =  <X(+)> 


Jk 


[nx  +  iny  +  %(p20  *  iZun  -  u02)] 

0  , 


<Y(+)>  and 

(»-nx)j(Y-ny)k) 


In  our  case,  X  =  0  and  H(+)  =  2tt(W(+)  -  W(-))  and: 

\  ■  (*<>))  -  o 

m20  =  0 


MU  = 


<eilM>  ■  (e 


i2tt  (Wd  (+)  -  W  (-)  )\  in„  -  tSU 


J>= 


Rv  '  "Rk  "\  _  e^"y  k02 


For  our  problem. 


(2n(W(+)  -  »(-)))  =  <2irft(+))  -  (2ttWC-))  =  0, 


(53) 


(54) 


since  in  a  stationary  process  the  expected  value  of  W(+)  and  W(-)  are 
identical.  We  also  find: 


J02  =  ([2lT(V+)  -  V-))  -  2*V+)  -  2ttWr(-)]2) 


=  2tt2[<Wr2(+))  -<2  \(+)\(-))  +  (\2(-))] 


2R(C,n') 
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R(£',n'))] 


(55) 


=  2tt2[2  (aw2  + 

=  2  (2tt)  2  cjw2[1  -  C  Cl '  ' )  ] , 


where 


cci'.n') 


- * — 2 - 

rt 


Here  C(5',n')  is  the  normalized  autocorrelation  function.  Equation  (52)  can 
then  be  written: 


<■ 


i2-rr(WR(+)  - 


V- 


!) 


~ (2tt) 2  5W2[1  -  C(£',n’)] 


where 


-  (2tt)2Sw2 


C’ 


(2^)25w2c 

~  (2tt)  2c^2 
e  -  1 


(57) 

(58) 


(59) 


When  the  irregularities  have  a  Gaussian  distribution  as  we  have  assumed,  we 
can  then  say  several  things  with  respect  to  Eq.  (58). 

(a)  We  have  placed  no  restrictions  (yet)  on  the  magnitude  of  o^. 

(b)  The  effective  transfer  function  can  be  thought  of  in  two  parts;  a 
scattered  part  with  an  arbitrary  autocorrelation  function  C',  and  a 
constant  specular  component.  These  components  are  illustrated  in 
Fig.  5a. 

(c)  The  associated  spread  function  is  a  delta  function  with  an  associ¬ 
ated  scattering  function.  The  Fourier  transform  of  C'  is  shown  in 
Fig.  5b. 

From  Eq.  (58)  we  see  that  the  total  flux  in  the  spread  function  distribu¬ 
tion  can  be  divided  into  two  parts,  that  is, 

$  =  $c  +  $H-  (60) 


18 


Fig.  5a.  Transfer  function.  Fig-  5b.  Spread  function. 


The  fraction  of  the  light  in  the  specular  part  is  given  by 


V* 


-(2^)25w2 


(61) 


and  the  scattered  fraction  by 


$  /$ 
H' 


1 


-(2tt)25w2 


(62) 


The  fractions  and  $(-.  can  be  determined  radiometrically,  and  the  magnitude 
of  the  irregularities  can  be  inferred  by  the  consequent  relation 


1 

(2r)2 


£n(l  +  $H/$C), 


which,  as  'S’h/ beconles  small,  approaches 


1 

C 2tt)  2 


c  w  • 


(63) 


(64) 
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The  shape  of  the  autocorrelation  function  of  the  irregularities  can  be 
inferred  from  the  measurable  scattering  function.  The  function  C'  in  Eq. 
(59)  and  the  scattering  function  are  a  Fourier  transform  pair,  and  the 
desired  autocorrelation  function  C  can  be  obtained  from  C'  by 

£n(l  +  (*/*)(]•) 

C  ■  SaTlyiffl  C6S) 

which,  for  small  becomes  simply 

C  x  C'.  (66) 

This  development  has  shown  that  it  is  possible  to  obtain  the  statistical 
properties  of  mirror  surface  irregularities  by  radiometric  measurement  of 
the  light  scattered  by  the  mirror. 
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Ill 


EXPERIMENTAL  APPARATUS  AND  TECHNIQUE 


Apparatus 

A'  schematic  diagram  of  the  apparatus  used  is  shown  in  Fig.  6.  To  mini¬ 
mize  the  effects  of  dust  and  dirt  the  equipment  was  set  up  in  a  closed-off 
area  of  the  Optical  Testing  Laboratory.  The  sample  and  the  associated 
detector  optics  were  further  isolated  by  means  of  a  black  velvet  tent  hung 
from  the  ceiling  and  enclosing  the  immediate  area. 

The  light  source  employed  was  a  Spectra-Physics  Model  165  Argon  Ion 
Laser.  The  laser  was  used  in  its  single-mode  configuration  with  the  power 
adjusted  to  20  mw.  Two  modes  of  regulation  are  available  in  this  laser; 
current  regulation  and  light  regulation.  In  all  instances  the  light  regula¬ 
tion  mode  was  used.  This  mode  varies  the  current  supplied  by  the  power 
supply  to  assure  intensity  regulation  to  within  one  percent.  A  further  fea¬ 
ture  of  this  particular  laser  is  an  aperture  within  the  cavity  which  permits 
adjustment  of  the  beam  diameter.  This  aperture  was  set  near  its  maximum 
opening,  300  pm,  to  minimize  the  effects  of  diffraction. 

Light  leaving  the  laser  was  directed  toward  the  sample  by  employing  two 
front-surfaced  mirrors.  The  beam  then  passed  through  a  chopper  which  modu¬ 
lated  the  beam  at  666  Hz.  The  chopped  light  then  passed  through  a  single 
blackened  aperture  before  striking  the  sample  surface.  Stray  light  associ¬ 
ated  with  the  laser  and  light  scattered  by  the  beam-directing  mirrors  and 
chopper  is  restricted  by  this  aperture.  The  sample  area  is  defined  by  the 
black  velvet  tent. 

The  specular  component  of  the  radiation  reflected  by  the  sample  passes 
back  through  the  tent.  In  the  case  of  "normal"  incidence  the  specular  beam 
returned  through  the  aperture  by  which  it  entered.  The  beam  in  such 
instances  is  not  precisely  normal  to  the  sample  surface  since  it  would  then 
effectively  become  part  of  the  laser  cavity  and  alter  the  characteristics  of 
the  incident  beam.  Instead  the  specular  and  reflected  beam  are  separated  by 
10  minutes  of  arc.  In  the  case  of  non-normal  incidence  the  specular  beam 
passes  through  a  well-defined  aperture  in  the  tent  and  is  subsequently 
directed  into  a  Rayleigh's  horn  light  trap. 

The  detector  in  all  instances  was  a  Phillip's  one- inch  end  on  photomulti¬ 
plier  having  an  S-20  photocathode.  Light  reaches  the  photomultiplier  by  way 
of  a  rigid  fiber-optic  probe.  Such  a  probe  offers  several  distinct  advan¬ 
tages  in  light  sampling.  First,  substantially  the  entire  cross-section  of 
the  fiber  bundle  is  a  collector.  In  addition,  the  field  of  view  can  be 
restricted  to  several  degrees  owing  to  the  propagation  properties  of  fiber 
bundles,  (see  Appendix  A),  and  lastly  the  diameter  of  the  bundle  was  small 
compared  to  that  of  the  photomultiplier  allowing  both  increased  angular 
resolution  throughout  the  scattered  field  and  the  ability  to  probe  within 
half  a  degree  of  the  incident  or  specular  beams. 
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Fig.  6.  Schematic  diagram  of  scatter  measurement  apparatus. 
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The  detector  probe  unit  was  mounted  on  a  rigid  arm  that  could  be  rotated 
in  either  of  two  orthogonal  directions.  Rotation  about  a  vertical  axis  was 
accomplished  by  means  of  a  massive  precision  rotary  table.  In  turn  the 
rigid  arm  was  attached  to  the  rotary  table  by  means  of  a  worm  gear  arrange¬ 
ment  which  allowed  rotation  about  the  horizontal  axis.  The  intersection  of 
these  two  axes  defined  the  sample  location.  Figure  7  illustrates  this 
configuration  in  detail. 

The  sample  holder  also  had  two  degrees  of  rotation;  one  about  the  verti¬ 
cal  axis  and  a  second  about  the  horizontal  axis.  The  horizontal  axis  always 
was  in  the  plane  of  the  sample  and  the  vertical  axis  always  intersected  the 
center  of  the  sample.  Alignment  of  all  four  axes  involved  in  sample  and 
detector  orientation  was  conducted  prior  to  each  experimental  run. 


Experimental  Procedures 


The  geometry  of  the  measurements  is  illustrated  in  Fig.  8.  The  distance 
from  the  sample  to  the  fiber-optic  probe  is  25  cm.  The  probe  is  .3  cm  in 
diameter  giving  it  an  angular  subtense  of  .25  degrees. 

In  order  to  facilitate  a  straightforward  data  reduction  routine,  only  the 
rotary  table  was  moved  during  any  data  run. 

The  measurement  procedure  is  then  as  follows:  the  laser  beam  is 
accurately  adjusted  to  strike  the  center  portion  of  the  sample  under  test  by 
moving  the  two  small  mirrors  previously  mentioned.  Subsequently  the  fiber- 
optics  probe,  sample,  and  rotary  table  are  properly  aligned  and  adjusted 
with  respect  to  the  incident  beam.  A  freshly  coated  magnesium  oxide  (MgO) 
reference  sample  known  to  have  a  diffuse  reflectance  in  excess  of  ninety- 
eight  percent  across  the  visible  spectrum  was  placed  in  the  sample  holder. 
Scattering  measurements  from  the  MgO  sample  were  integrated  to  calculate  the 
incident  intensity  of  the  incoming  radiation. 

The  sample  to  be  measured  is  then  placed  in  the  holder  and  the  scatter 
field  is  measured  by  moving  the  fiber  probe  from  a  position  near  the  specu¬ 
lar  beam  toward  the  plane  of  the  sample.  This  motion  is  accomplished  by 
moving  the  rotary  table  in  steps  toward  increasingly  larger  angles. 

Typically  the  scatter  was  measured  at  thirty  different  angles  during  a  given 
run.  At  each  angle  a  voltage  measurement  proportional  to  the  intensity 
level  of  the  scattered  field  was  recorded.  These  measurements  were  then 
recorded  as  the  ratio  of  scattered  light  to  reflected  light  vs.  angle  from 
the  specular  beam.  This  comprises  the  raw  data. 

To  compensate  for  unwanted  stray  light  the  sample  was  entirely  removed 
from  its  holder  and  the  beam  was  allowed  to  exit  through  the  black  velvet 
tent.  In  this  sense  the  sample  can  be  thought  of  as  perfectly  reflecting 
and  any  resulting  measurements  of  scattered  radiation  can  then  be  attributed 
to  background.  A  background  measurement  was  taken  in  conjunction  with  each 
run.  These  readings  were  subtracted  from  the  raw  data. 

Another  source  of  unwanted  light  was  due  to  the  divergence  and  finite 
extent  of  the  laser  beam  itself. 
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Compensation  for  the  error  introduced  by  beam  spread  was  made  in  a  simi¬ 
lar  fashion  to  that  used  for  the  background  measurements.  Again  the  sample 
was  removed  from  the  sample  holder  and  "no  sample"  was  considered  to  be  the 
perfect  reflector.  In  this  instance,  however,  the  probe  was  moved  to  the 
position  of  its  virtual  image  as  seen  when  looking  into  the  actual  sample. 
This  geometry  is  shown  in  Fig.  9.  The  probe  is  then  moved  away  from  the 
transmitted  beam  and  measurements  are  taken  at  angles  corresponding  to  those 
taken  when  the  actual  sample  was  in  place.  These  measurements  are  also  sub¬ 
tracted  from  the  initial  raw  data.  They  were  recorded  as  "beam-spread" 
measurements.  These  measurements  became  noise  limited  when  the  probe  was 
6  degrees  from  the  incident  beam. 


( 


3  j\  l . . 1 

/  |  Sample  Holder 
1 

/ 


C> 

Detector 


Incident  Beam 


i 


Fig.  9.  Geometry  of  the  beam  spread 
measurements. 
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IV 


DATA  AND  RESULTS 


The  distribution  of  light  scattered  by  a  plane  mirror  as  a  function  of 
surface  roughness,  angle  of  incidence,  and  wavelength  are  reported  in  this 
section.  General  results  are  compared  with  the  theory. 


Background  Data 


As  noted  in  the  section  on  Experimental  Procedure,  the  background  mea¬ 
surements  were  due  to  spread  of  the  laser  beam.  As  these  measurements  were 
all  taken  through  the  sample  holder  with  no  sample  in  place,  it  was  expected 
that  all  measurements  would  be  the  same.  Figures  10  and  11  show,  however, 
that  the  values  of  the  background  varied  considerably.  The  cause  for  these 
background  differences  is  not  clear,  but  several  things  should  be  pointed 
out  with  respect  to  them.  First,  the  angular  extent  of  the  background  is 
confined  to  within  eight  degrees  of  the  specular  beam,  with  most  detail  of 
any  note  being  within  six  degrees  of  the  specular  beam.  In  some  cases  the 
small  angle  data  has  a  shape  very  similar  to  that  of  a  particular  background 
measurement,  and  in  others  there  is  no  strong  correlation.  The  net  result 
is  that  these  background  measurements,  which  should  all  be  identical,  reduce 
our  confidence  in  the  small  angle  data  (that  data  between  1.5°  and  6°). 

Data  taken  outside  this  range  appears  to  be  reliable. 

Because  of  discrepancies  in  the  background  measurements  we  have  presented 
much  of  the  data  in  a  dual  format.  In  the  first  format  we  subtracted  the 
background  that  was  taken  in  conjunction  with  the  particular  set  of  measure¬ 
ments.  The  second  format  presents  the  data  with  the  "best-fit"  background 
subtracted . 


Presentation  of  the  Data 


A  measure  of  the  total  flux,  3>,  incident  on  the  samples  was  obtained  by 
performing  a  numerical  integration  of  the  scattered  light  distribution  from 
the  magnesium  oxide  reference  sample.  These  integrations  were  made  in  con¬ 
junction  with  each  set  of  measurements. 

Background  data  describing  the  flux  distribution  of  the  incident  beam  is 
then  subtracted  from  measurements  of  the  intensity  of  the  scattered  field. 
These  data  are  represented  by  n ' .  The  corrected  data  is  then  normalized  by 
dividing  by  $  and  the  data  from  each  run  can  now  be  plotted  to  yield  the 
scattered  light  distribution  of  the  sample. 


These  distributions  are  plotted  as  modified  log-log  plots  with  the 
ordinate  being  log  n/cos9  (where  n  =  o'/^  is  the  scatter  coefficient);  where 
0  is  the  angular  position  of  the  measurement  with  respect  to  the  normal  to 
the  sample.  The  abscissa  in  these  plots  is  log  (sin9  -  sin  i) ;  where  i  is 
the  angle  of  incidence  (see  Fig.  8) .  We  now  have  0  =  k+i,  k  being  the  angle 
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from  the  specular  beam  to  the  measurement.  To  illustrate  the  effect  that 
these  plotting  coordinates  exercise  on  the  data,  intermediate  steps  arc 
illustrated  in  Figs.  12,  13,  and  14  for  samples  247  and  MgO. 

Figure  12  is  a  log-log  plot  of  n  vs.  the  angle  k  from  the  specular 
beam.  In  Fig.  13  we  have  plotted  the  same  information,  however,  we  have  now 
divided  the  scattered  flux  by  cos0.  The  effects  of  this  division  are  seen 
easily  at  the  larger  angles  where  cos0  is  small  and  the  plotted  values  are 
noticeably  higher.  The  result  of  this  division  is  of  particular  note  in  the 
case  of  the  MgO  sample  where  the  curve  has  become  almost  flat.  (A  flat  plot 
would  indicate  a  perfect  Lambertian  surface.) 

Figure  14  illustrates  the  change  of  scale  in  the  abscissa  when  we  shift 
from  log  (k)  to  log  (sin6  -  sin  i) .  In  this  case  we  have  used  data  for 
i  =  0°  and  i  =  45°  from  sample  247.  The  change  in  scale  for  normal  inci¬ 
dence  measurements  can  be  noted  for  both  MgO  and  sample  247. 


The  data,  at  least  in  the  large-angle  region  where  it  is  well-behaved, 
appears  to  follow  an  inverse  power  law,  and  must  eventually  level  off  to  a 
finite  maximum  value  for  very  small  angles,  although  this  leveling  does  not 
show  in  the  data.  A  function  that  may  be  used  to  represent  this  behavior  is 


i  =  io[i  +  (p/po)2r(s/2). 


(67) 


where  I0  is  the  central  maximum  radiant  intensity,  p  =  /a2+32  is  the  direc¬ 
tion  cosine  of  the  scattering  angle  (assuming  the  scattering  function  is  a 
figure  of  revolution) ,  pD  is  a  scale  factor  measuring  the  width  of  the 
scattering  function,  and  s  is  the  rate  at  which  the  function  decreases  for 
large  angles.  For  small  angles  where  p/pq  <<  1> 


1  -  V1  -  f  Cp/p0)2]> 

and  for  large  angles  where  p/p0  >>  1, 

I  =  (Ip  S)p'S. 
oo 

The  total  flux  scattered  out  to  a  radius  p*  is  given  by 

P* 

<!>s(f>*)  =  2rr  j  J  p dp 

0 

and  the  total  flux  scattered,  $  (1)  =  is  obtained  by  setting  p*  = 


2  it  I  p  2 

o  o 

2-s 


[1  +  (P*/P0)2J 


2,%(-S) 


!> 


(68) 


(69) 


s  <  2, 
(70) 
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In  the  limit  for  pQ  «  l,  we  obtain 


2  ttI  p  S 
oo 

H  ~  2- 


(71) 


If  we  compare  this  with  Eq.  (69),  we  see  that  we  can  obtain  an  approximate 
value  for  the  total  scattered  flux  by  the  relation 


$ 


H 


2tt 

2^s 


Id), 


(72) 


without  knowing  the  values  of  I  and  p  . 

o  o 


Scatter  as  a  Function  of  Surface  Roughness 


Variations  in  surface  roughness  were  obtained  by  employing  a  variety  of 
grinding  polishing  techniques  as  outlined  in  Appendix  A.  The  scattered 
light  distributions  of  all  samples  listed  in  Table  3,  Appendix  A,  were  mea¬ 
sured  at  normal  incidence.  The  results  of  these  measurements  are  shown  in 
Figs.  15  and  16. 

It  should  be  noted  that  the  slope  of  the  wide  angle  portion  of  the  pol¬ 
ished  samples  is  close  to  being  flat  for  the  angles  beyond  six  degrees.  The 
slopes  vary  from  about  1.25  to  2.0.  Also  of  significance  is  the  change  of 
shape  brought  about  by  subtracting  the  appropriate  background.  Selection  of 
this  background  appears  to  have  substantially  removed  the  anomaly  in  the 
plots  of  samples  236,  222,  247,  and  EOF  3-4. 

Samples  221  and  229  are  ground  samples  and  exhibit  a  behavior  close  to 
that  of  MgO.  No  specular  beam  was  detectable  when  these  samples  were  illu¬ 
minated  at  normal  incidence. 

Separate  scattering  phenomena  may  account  for  the  kink  in  the  plots  of 
polished  samples  at  approximately  six  degrees.  However,  uncertainty  of  the 
background  data  necessitates  further  research  to  determine  if  this  is  true. 


Scatter  as  a  Function  of  Angle  of  Incidence 


In  the  general  case,  light  is  not  normally  incident  upon  an  optical  sur¬ 
face.  The  behavior  of  scattered  radiation  when  light  is  incident  upon  a 
sample  at  angles  between  0°  and  60°  was  investigated.  In  all  cases  the  mea¬ 
surements  were  made  in  the  plane  of  incidence  and  only  include  measurements 
in  the  shaded  portion  of  Fig.  8.  (A  modification  to  the  existing  apparatus 
is  suggested,  [see  "Suggestions  for  Future  Research"]  that  will  allow  mea¬ 
surement  throughout  the  plane  of  incidence  and  planes  orthogonal  to  the 
plane  of  incidence.) 
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DISTANCE  FROM  SPECULAR  BEAM  (sin  0  -  sin  i) 


Fig.  16.  Scatter  as  a  function  of  surface  roughness,  normal  incidence; 
corrected  for  best  fit  background  (see  discussion) . 
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The  results  of  a  series  of  measurements,  with  the  beam  incident  at  vari¬ 
ous  angles,  is  shown  in  Figs.  17  and  18-  In  Fig.  17  a  strong  anomaly  is 
present,  which  is  substantially  removed  when  a  more  judicious  background  is 
selected  as  shown  in  Fig.  18. 

The  slopes  of  these  curves  is  close  to  1.5.  In  Fig.  18  the  individual 
curves  exhibit  close  agreement  with  each  other.  More  extensive  research 
should  be  pursued  to  substantiate  this  type  of  agreement. 


Scatter  as  a  Function  of  Wavelength 

An  original  intent  of  the  research  reported  here  was  to  measure  the 
scattered  light  distribution  of  samples  at  two  different  wavelengths:  one 
in  the  visible  portion  of  the  spectrum  and  a  second  in  the  mid-infrared.  As 
research  progressed,  it  became  evident  that  satisfactory  results  could  not 
be  obtained  if  the  original  wavelength  regions  were  retained.  A  decision 
was  made  to  use  two  wavelengths  in  the  visible  portion  of  the  spectrum  for 
the  following  reasons: 

(1)  Two  measurements  in  such  widely  separated  regions  of  the  spectrum 
would  require  that  the  separate  laser  sources  be  used.  It  is 
unlikely  that  these  sources  could  be  aligned  to  illuminate  identical 
portions  of  the  sample  under  test. 

(2)  In  order  to  measure  the  light  levels  associated  with  wide  angle 
scatter,  two  separate  detectors  are  necessary.  It  is  difficult  to 
insure  that  the  absolute  energy  and  the  collection  area  of  each 
detector  be  well  known. 

(3)  A  reference  sample  that  would  be  useful  in  both  the  infrared  and 
visible  was  unknown. 

In  contrast,  the  availability  of  an  argon-ion  laser  simplified  the  associ¬ 
ated  problems  of  measurement  immediately.  A  prism  internal  to  the  cavity  of 
the  laser  permits  selection  of  several  laser  lines  without  the  alteration  of 
beam  direction.  The  lines  used  were  514.5  nm  (green)  and  459.0  nm  (blue). 

A  photomultiplier  exhibits  good  response  at  each  of  the  available  lines 
allowing  a  single  detector  configuration  to  be  used  for  both  wavelengths. 
Furthermore,  the  diffuse  reflectance  of  MgO  exceeds  .98  throughout  the  spec¬ 
tral  region  of  the  laser.  This  allowed  an  MgO  sample  to  be  used  as  a  refer¬ 
ence.  Detector  response,  laser  power  output,  reflectance  of  intermediary 
mirrors,  and  transmittance  of  the  fiber-optic  probe  could  all  be  reduced  to 
a  single  proportionality  factor.  This  factor  is  obtained  by  taking  the 
ratio  of  the  integrated  diffuse  reflectance  of  the  MgO  reference  at  the  two 
wavelengths  of  concern. 

To  measure  relative  scatter  as  a  function  of  wavelength,  an  MgO  reference 
sample  is  placed  in  the  sample  holder.  Its  diffuse  reflectance  is  measured 
at  the  first  of  the  two  wavelengths.  The  reference  is  removed  from  the 
sample  holder  and  replaced  with  the  sample  under  study.  The  diffuse  reflec¬ 
tance  of  the  sample  is  then  measured  at  the  first  wavelength.  The  laser  is 
then  tuned  to  the  second  of  the  chosen  wavelengths,  and  the  diffuse 
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reflectance  is  measured  at  this  wavelength.  (It  should  be  noted  that  the 
change  of  wavelength  does  not  alter  the  position  of  any  optics  outside  the 
laser  cavity,  and  that  the  laser  beam  in  both  cases  strikes  the  same  portion 
of  the  sample.)  The  test  sample  is  now  replaced  by  the  reference  sample  in 
the  sample  holder,  and  the  diffuse  reflectance  of  the  reference  at  the 
second  wavelength  is  measured. 

A  numerical  integration  performed  on  the  data  taken  from  the  reference 
measurements  yields  the  relative  intensity  of  the  incident  radiation  with 
the  thru-put  of  the  entire  system  (reflectance,  transmission,  photocathode 
response  and  electronics)  taken  into  account.  The  ratio  of  the  resultant 
integrations  is  a  scale  factor  that  is  then  used  to  "correct"  the  measure¬ 
ments  from  the  test  sample  for  all  effects  except  for  wavelength.  The  re¬ 
sult  of  these  measurements  is  shown  in  Figs.  19,  20, $21.  All  measurements 
were  made  at  normal  incidence. 

In  order  to  predict  the  wavelength  dependence  of  the  scattering  function 
for  large  angles,  we  first  combine  Eqs.  (69)  and  (71),  obtaining 

1  =  ™ 

The  form  of  the  scattering  function  should  not  vary  with  wavelength,  so  the 
slope  s  should  not  be  dependent  on  wavelength.  Furthermore,  if  the  scat¬ 
tered  fraction  of  the  flux  is  small,  the  total  scattered  flux  is  given  by 

-(2tt)25  2 

$H  =  1  -  e  W  ~~  (2tt )2(ow/X)2,  (74) 


so  that  $  ,  and  therefore  I,  should  vary  inversely  as  the  square  of  the 
wavelength.  This  is  consistent  with  the  data. 


Numerical  Integration  of  the  Scattered  Light  Distributions 


A  numerical  integration  of  the  scattered  light  distribution  yields  4>h, 
the  net  scattered  flux  if  we  can  assume  the  distribution  is  rotationally 
symmetric,  and  that  the  contribution  due  to  the  fraction  of  light  not 
included  in  the  integral  is  small.  These  assumptions  appear  reasonable  as 
we  have  no  evidence  to  the  contrary.  The  results  of  these  integrations  are 
listed  in  Table  1. 

Further,  from  the  theory  we  know  that  the  rms  surface  roughness,  a  . 
is  given  by 
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Fig.  19.  Scatter  as  a  function  of  wavelength,  sample  #236. 
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Fig.  20.  Scatter  as  a  function  of  wavelength,  sample  #247. 
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Since  $  =  %+$(>  we  can  ma^e  the  assumption  that  =  4>  when  3>^  is  very 

small.  This  is  the  case  for  all  but  the  two  ground  surfaces,  samples  221 
and  229.  The  rms  roughness  is  then  given  by 


(76) 


Calculations  of  ow  are  given  in  Table  2.  These  results  exhibit  an  order  of 
magnitude  agreement  with  the  estimated  roughness  values. 


TABLE  1 

RESULTS  OF  THE  NUMERICAL  INTEGRATION 
OF  THE  SCATTERED  LIGHT  DISTRIBUTIONS 


Sample  No.  Integrated  Scatter  (%) 


236 

0.014 

247 

0.033 

EDF 

3-4 

0.038 

222 

0.046 

EDF 

3-11 

0.38 

229 

69.5* 

221 

71.1* 

*In  the  case  of  the  ground  samples,  no 
specular  beam  was  apparent  and  these 
figures  represent  the  total  diffuse 
reflectance . 


TABLE  2 


CALCULATED  RMS 

SURFACE  ROUGHNESS,  a 

w 

Sample  No. 

$H/<i> 

°w  (nm) 

236 

0.00014 

0.95 

247 

0.00033 

1.5 

EDF  3-4 

0.00038 

1.6 

222 

0.00046 

1.8 

EDF  3-4 

0.0038 

5.1 
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SUGGESTIONS  FOR  FUTURE  RESEARCH 


In  conducting  the  research  presented  in  this  report  several  ideas  for 
improving  the  measurements  and  our  understanding  of  surface  scattering  have 
been  kindled.  These  ideas  are  summarized  below. 

(1)  A  redesign  of  the  existing  apparatus  is  possible  that  will  allow  one 
to  scan  through  the  specular  beam  and  thus  make  measurements  over 
the  entire  plane  of  incidence.  This  can  be  accomplished  by  intro¬ 
ducing  an  offset  into  the  fiber-optic  probe. 

(2)  Through  the  use  of  a  lens  to  focus  the  incident  laser  beam  onto  the 
hemisphere  an  exact  Fourier  transform  can  be  formed.  Approximations 
involving  the  quadratic  phase  factor  can  then  be  eliminated. 

(3)  An  investigation  of  surface  profiles  using  the  techniques  of  stereo¬ 
electron  microscopy  could  yield  more  accurate  measurement  of  auto¬ 
correlation  function  and  the  surface  height  distribution. 

(4)  Surfaces  that  are  partially  polished,  and  surfaces  that  are  polished 
and  then  partially  ground  as  well  as  materials  other  than  glass 
should  be  investigated. 


'  ) 
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APPENDIX  A 


SAMPLE  PREPARATION 


This  section  deals  with  sample  selection,  polishing,  cleaning  and 
coating.  An  original  goal  of  the  research  reported  here  was  to  choose  sam¬ 
ples  that  bracketed  the  types  of  surfaces  typically  produced  in  optical 
component  manufacture.  Coatings  for  these  surfaces  we  produced  were  care¬ 
fully  controlled  and  in  all  cases  we  believe  had  less  effect  on  overall 
scatter  than  did  the  condition  of  the  substrates.  All  surface  roughnesses 
are  reported  as  a  combination  of  coating  plus  substrate  roughness. 


Substrate  Materials  and  Preparations 

Samples  used  in  these  experiments  were  of  two  different  types  of  glass. 
The  hardness  range  of  conventional  optical  glass  was  well  bracketed  by  the 
choice  of  substrates.  The  large  majority  of  the  substrates  were  made  of 
fused  silica  of  the  Homosil  trademark.  These  substrates  were  of  exceptional 
quality  and  essentially  free  from  bubbles  and  sleeks.  All  substrates  were 
3.17  cm  in  diameter  and  0.63-cm  thick.  Fused  silica  is  one  of  the  hardest 
of  the  conventional  optical  substrates.  The  remaining  substrates  were  made 
of  optical  quality  extra  dense  flint  glass.  The  particular  category  in  this 
case  is  known  as  EDF3.  EDF3  is  one  of  the  softest  optical  glasses  and  it 
was  chosen  as  a  substrate  for  this  particular  reason. 

Surface  preparation  of  the  different  samples  is  outlined  in  Table  3.  All 
samples  were  finished  to  be  closely  flat.  Grinding  and  polishing  procedures 
are  more  closely  described  by  Mott  (1971)  and  Orme  (1972)  .  Surface  rough¬ 
ness  in  the  samples  varied  from  1.0  nm  to  3,000  nm.  The  EDF3  samples  were 
much  more  prone  to  small  scratches  as  might  be  expected.  In  all,  thirty-two 
samples  were  prepared  and  of  these,  eight  were  selected  for  investigation. 


Cleaning 

The  prepared  samples  were  cleaned  prior  to  coating  with  aluminum. 

Cleaning  consisted  of  careful  washing  with  Liquinox,  a  mild  detergent,  under 
very  warm,  filtered  tap  water.  Such  soaping  was  repeated  three  times  for 
each  individual  sample.  Samples  were  then  mounted  in  a  sample  holder  while 
held  in  distilled  water.  Once. in  a  holder,  samples  were  moved  to  an  ultra¬ 
sonic  cleaner  filled  with  distilled  water  for  rinsing.  Purity  of  the  dis¬ 
tilled  water  in  the  ultrasonic  cleaner  was  monitored  through  resistance  mea¬ 
surements.  The  water  was  changed  several  times,  until  it  was  substantially 
as  pure  as  it  was  before  the  samples  were  added.  Once  rinsed  the  samples 
were  set  to  dry  in  a  dust-free  atmosphere.  Dry  samples  were  removed  from 
the  sample  holders  and  placed  in  individual  boxes  being  supported  by  the 
edge  of  their  backside.  Mott  (1971)  used  a  similar  cleaning  technique, 
which  he  describes  more  completely. 


Coating 


Cleaned  dry  samples  were  then  placed  in  a  high  vacuum  chamber  and  coated 
to  opacity  with  pure  aluminum.  Coating  technique  varied  from  the  standard 
only  in  that  excessive  care  was  taken  to  allow  the  chamber  to  reach  a  pres¬ 
sure  below  2  x  10"r>  torr  prior  to  coating.  Opacity  was  monitored  by 
watching  the  extinction  of  the  tungsten  evaporating  filament  through  the 
substrates.  The  samples  were  allowed  to  cool  to  room  temperature  prior  to 
removal  from  the  chamber.  Each  coating  run  contained  ten  different  samples. 
Once  coated,  samples  were  returned  to  their  individual  storage  boxes.  After 
all  samples  were  coated,  the  best  samples  of  each  type  were  selected  for 
measurement.  This  selection  was  made  on  the  basis  of  individual  inspection 
of  each  sample  while  held  under  a  microscope  illuminator  in  an  otherwise 
dark  room.  Samples  that  had  coating  nonuniformities,  sleeks  or  pinholes 
were  rejected  as  were  those  with  waterspots,  large  scratches  of  otherwise 
questionable  appearance.  Prior  to  each  scatter  measurement  run,  samples 
were  again  individually  inspected  for  flaws.  Dust  was  removed  using  a 
commercially  available  pressurized  air  can.  After  each  sample  was  measured 
for  scatter,  it  was  returned  to  its  individual  box. 
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TABLE  3 


SURFACING  PROCEDURES  AND  ROUGHNESSES 


Sample 

Number 

Substrate 

Material 

Lap 

Material 

Polishing 

Medium 

Polishing 

Time 

Force  Applied 
kg/ cm2 

Approximate 
Surface  (nm) 
Roughness 

236 

Quartz 

Optical  Quality 
Pitch 

CeO*(milled  2000  hr) 

8  hr 

.2 

1.0 

Distilled  H^O*  w/same 

lap 

8  hr 

.2 

247 

Quartz 

Wood  (end  cut 
fir) 

CeO  (unmilled) 

4  hr 

.3 

4.0 

222 

Quartz 

Plate  Glass 

A^O^Clym  diam) 

1  hr 

.12 

12 

221 

Quartz 

Cast  Iron 

A^O^^nm  diam) 

20  min 

.15 

100 

229 

Quartz 

Cast  Iron 

A^O^^Oym  diam) 

20  min 

.  15 

1500 

EDF3-4 

EDF-3 

Optical  Quality 
Pitch 

CeO (mi lied  2000  hr) 

8  hr 

.  12 

2.5** 

H^O* (distilled  w/same 

lap) 

8  hr 

.  12 

EDF3-11 

EDF-3 

Optical  Quality 

CeO  (milled  500  hr) 

4  hr 

.12 

15** 

Pitch 


^Closed  circulation  system. 

**Soft  glasses  exhibit  a  strong  tendency  to  scratch  and  sleek. 


APPENDIX  B 


PROPERTIES  OF  THE  FIBER-OPTIC  PROBE 


The  fiber-optic  bundle  offers  several  unique  advantages  in  the  collection 
of  the  scattered  radiation.  Being  rigid,  it  is  self-supporting  and,  in 
effect,  transfers  the  detector  to  the  end  of  the  fiber  bundle.  Furthermore, 
being  small  in  cross-section,  one  is  allowed  to  probe  the  space  near  the 
specular  direction  without  the  detector  scattering  light  back  towards  the 
sample.  Finally,  fiber-optic  bundles  have  the  property  that  off-axis  light 
leaves  the  bundle  both  attenuated  and  distributed  away  from  the  normal 
direction.  Figures  22  and  23  illustrate  this.  Correct  baffling  between  the 
end  of  the  fiber-optic  bundle  and  the  detector  can  then  eliminate  much  of 
the  stray  and  background  radiation. 
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Fig.  22.  Relative  light  output  vs.  input  direction  of  collimated 
light  for  a  fiber-optic  probe. 


OUTPUT 


Fig.  23.  Relative  output  intensity  vs.  output  direction  for  a 
fiber-optic  probe. 
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